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Abstract 

Simultaneous measurements of 
surface atmospheric parameters 
and cloud liquid water are used 
to test and compare the 
accuracy of three different 
cloud models. 


1. Brief Review of Cloud 
Attenuation Models 

Numerous models for predicting 

the attenuation of electro- 
magnetic waves propagating 
through clouds were developed 
over the years from a variety 
of theoretical and empirical 
methods. Cloud modeling for 
the purposes of assessing 
attenuation can be divided into 
essentially three different 
catagories: 1) attenuation is 
computed by using a Rayleigh 
approximation to Mie scattering 
theory [Gunn, East, 1954], 
[Staelin 1966], [Liebe, Manabe, 
Hufford, 1989]; 2) attenu- 
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ation is directly correlated to 
surface absolute humidity 
[Altshuler, Marr, 1989]; 3) 
Meteorological data and compu- 
tations are used to determine 
cloud liquid and then attenu- 
ation is computed using a 
slightly modified version of 
the catagory 1 models described 
above. [Slobin, 1982], 
[Dintelmann, Ortgies, 1989] . 

Although their math- 
ematical form and predictions 
vary over a fairly large range, 
a parameter common to all 
models is the liquid water 
content of the cloud. Unfor- 
tunately, this fundamental 
parameter is also the most 
difficult to predict and to 
measure. 

A detailed comparison of 
five prominent cloud models 
developed over the last forty 
years shows good agreement at 
frequencies below 40 GHz for 
light to medium clouds cond- 
itions [Gerace, Smith, 1990] . 
However, for heavy to very 
heavy clouds and frequencies 
above 10 GHz, the models 
diverge from each other. 

The recent availability of 
radiometric measurements of 
atmospheric parameters and the 
worldwide availability of 
surface atmospheric meas- 
urements have inspired the 
development of new cloud at- 
tenuation models. These new 
models strive to relate surface 
atmospheric measurements to 
cloud attenuation. The overall 
underlying assumption is that 
the liquid water content of 
clouds is in some way related 
to the water vapor present at 
the earth's surface. 

This paper describes how 
three of these new cloud models 
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perform on a cloud event that 
was in no way related to the 
empirical data used to develop 
the models. The preliminary 
results presented below are an 
attempt to qualitatively verify 
both the mathematical cloud 
models (types 2 and 3) and the 
latest methods available for 
extracting data from an 
independent cloud event. A 
complete statistical analysis 
is forthcomming when we 
complete our analysis using 
cloud data measured au numerous 
sites worldwide. 

We begin by introducing 
the Altshuler-Marr, Dintelmann- 
Ortgies, and GSW cloud 
attenuation models and briefly 
discussing a method for 
measuring cloud liquid water. 
We then present our methods for 
comparing the models along with 
graphical results. The results 
are also cross checked with the 
well established Slobin cloud 
models [Slobin, 1982] . 

2 . Altshuler Model 

By correlating data of 
absolute surface humidity with 
measurements of zenith cloud 
attenuation in the Boston area, 
Altshuler derived the following 
empirical model for a nominal 
cloud temperature of 10'C 
[Altshuler, 1989] : 


«=|-0.0242+0. 000751 + -^^ {11.3+p) 

( 1 ) 

where 

a=zenith attenuation (dB) 
X=wavelength (mm) 
p=surface water vapor 
density (g/m 3 ) 


To account for elevation angles 
other than 90 degrees, eq. 1 
must be multiplied by the 
following: 

c\sc(0) 

^) =v (a e +h e ) 2 -a e 2 cos 2 (Q)] 2 
-a e sin(0) 

( 2 ) 

where 

0 * elevation angle 
a e = effective radius of 
the earth (4/3 earth taken as 
8497 km) 

h. = 6.35-0.302p effective 
cloud height (km) 

p - surface absolute 
humidity (g/m 3 ) 

While the Altshuler model 
is primarily an empirical model 
the next model is more appro- 
priately classified as a 
semiempirical model. 


3 . Dintelmann-Ortgies Model 

Using standard meteoro- 
logical equations along with 
radiometer attenuation and 
concurrent meteorological 
measur e ments, Dintelmann and 
Ortgies derived the following 
semiempirical model for cloud 
attenuation prediction 
[Dintelmann, Ortgies, 1989] : 

«= p o^(i-~ (g/m 3 ) 

(3) 


where 
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M = cloud liquid water 
(g/m 3 ) 

T 0 = surface temperature 

T = cloud temperature 

p 0 = surface water vapor 
density (g/m 3 ) 

K = the ratio of the 
specific heat of water at 
constant pressure to the 
specific heat of water at 
constant volume (approximately 
= to 4/3) 

g=acceleration of gravity 
(9.8 m/s 2 ) 

R=gas constant for air 
(approximately 287 J/K-Kg) 

H = height of the 0 degree 
isotherm (m) 

The height of the 0 degree 
isotherm can be approximated 
by : 


H- 0 . 89 +0 . 165(T 0 -273) 

(4) 


where 

T 0 = Surface Temp. (K) 

Then the attenuation 
through the cloud can be 
computed using an equation 
Dintelman borrowed from 
[Slobin, 1982]: 

4 . 343*10°- 0122 < 291 -r> -i 
« = — ■ J xu -1 .16Af 


(5) 

where a is now in dB/km, T is 
the cloud temperature in 
Kelvin, and X is the wavelength 
in centimeters. 

To obtain the total attenuation 
through the cloud, Dintelman 
used radiometer measurements to 
obtain the following empirical 
formula for the cloud vertical 


extent : 


A^.lS-O.OZBW+O.OOSSAf 2 (km) 

( 6 ) 


where M is the cloud liquid in 
g/m 3 . 

Inherent in this model is 
the assumption that clouds form 
around the O'C isotherm. The 
next model attempts to refine 
the Dintelmann-Ortgies model by 
including a calculation aimed 
at predicting more accurately 
the altitude of cloud 
formation. 

4. GSW Modal 

The altitude at which the 
actual water vapor density 
exceeds the saturated water 
vapor density for the 
temperature and pressure at 
that point is called the 
lifting condensation level. The 
GSW (initials of authors' last 
names) model assumes that this 
is the altitude at which clouds 
begin to form. The model can 
be described as follows: 

The initial version of the 
GSW model assumes a linear 
adiabatic temperature lapse 
rate of 6 deg C per kilometer: 

T(h) =T a -yT 

y=6° / Km 


(7) 


Then a vertical saturated 
water vapor profile can be 
computed as follows: 
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( 8 ) 



Now a vertical water vapor 
profile can be computed as 
follows : 


where e, is the water vapor 
pressure and is given by the 
following formula due to 
[Nordquist, 1973] : 


(13) 


e e =lO x where 


x^c^l. 3816e' 7 ' loPl 

+8.1328e' 31oP2 

_ 2949.076 
T 


( 9 ) 


where 

p 1 =ll. 344 -0. 030399 8T 

( 10 ) 


One can compute the 
lifting condensation level by 
equating equations (8) and (13) 
and solving for the height, h. 
This is where the saturation 
vapor density equals the actual 
vapor density and is most 
likely the altitude at which 
the cloud begins to form. 

Above the lifting 
condensation level, water vapor 
continues condensing as long as 
the actual vapor density 
exceeds the saturated vapor 
density. Loosely based on 
actual measurements of total 
integrated cloud liquid water 
and typical values of cloud 
liquid water densities, an 
estimate of the cloud liquid 
water content can be computed 
as follows: 


Af=p (A) -p s {A) 


Pa =3. 49149- 1302 r 8844 

(ID 


c x = 23 . 832241-5 . 02808logCT) 

( 12 ) 


( 14 ) 


where h' is the altitude at 
which p = 1.25 p, . 

Then cloud attenuation can 
be computed using equations 5 
and 6 with equation 6 modified 
by multiplying all of the 
coefficients by a factor of 
ten. This factor of ten will 
most likely be refined as we 
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average in more data sets from 
various sites to improve our 
model . 

Next, we describe a method 
for measuring the amount of 
liquid water in a cloud. 

5 . Cloud Liquid Hater 
Measurements 

Radiometer measurements of 
atmospheric absorption at two 
frequencies, a water vapor 
sensitive frequency and a cloud 
liquid water sensitive 
frequency (say 20.6 and 31.65 
MHz), can lead to a 
determination of total 
integrated cloud liquid water, 
L [Westwater, 1978] . The 
computation can be summarized 
as follows: 


£- ( K VU-^l' fK V2-^u) 

(K W *Lu-KvuK£l) 

( 15 ) 


K m = path averaged 
absorption coefficient of 
liquid at the upper liquid 
water sensitive frequency, u. 

K vi = path averaged 
absorption coefficient of vapor 
at the lower water vapor 
sensitive frequency, 1. 


K Li f * path averaged 
absorption coefficient of 
liquid at the lower liquid 
water sensitive frequency, 1. 

T mr = mean radiating 
temperature 

T b b = cosmic background 
"big bang" brightness 
temperature (2.8 K) 

T bv = measured value of 
the microwave brightness 
temperature at frequency, v. 

T dv = dry absorption at 
frequency, v. 


where 


( T mi~Tt ft) 

( 16 ) 


for v = 1, u 


where 

k Vu = path averaged 
absorption coefficient of vapor 
at the upper liquid water 
sensitive frequency, u. 


Measurements of cloud 
liquid water using the above 
algorithm are currently being 
made by the Wave Propagation 
Laboratory (WPL) of the 
National Oceanic and 
Atmospheric Administration 
(NOAA) at San Nicolas Isand, 
CA, and Denver Colorado. 

We are now intensively 
analyzing data that was 
collected throughout the 1980s. 
The results in this report are 
based on data taken in July 
1984. 

6. Method of Comparison 

Figure 1 depicts our 
method of comparison. Using 
surface atmospheric 
measurements taken in Denver 
CO, cloud liquid water contents 
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Figure 1. Method of Comparison 


computed using the Dintelmann- 
Ortgies and GSW models were 
compared to measurements of 
cloud liquid present at the 
time the surface measurements 
were recorded. Attenuation 
predicted by the Altshuler 
model was compared to that 
predicted by the Dintelmann- 
Ortgies model (via the Slobin 
approximation discussed above) . 


7 . Results . 

A time series of the 
surface measurements taken 
during a sample cloud event is 
shown in figure 2. Figure 3 
shows a comparison of the 
Dintelmann predictions to 
NOAA' s measurements of cloud 
liquid water. Figure 4 shows a 
similar comparison for the GSW 
model. Note that the order of 
magnitude of the total 
integrated liquid (cm) for all 


three models is correct. 
However, the shape of the 
curves agree qualitatively only 

during the last half of the 
three hour measurement period. 
Also note that the Dintelmann- 
Ortgies model predicts high 
liquid water content (g/m 3 ) and 
low vertical cloud extent as 
compared to the Slobin models 
described in figure 5. But the 
two effects sort of cancel each 
other out when computing the 
total integrated liquid (cm) 
because the units conversion 
from g/m 3 to cm is as follows: 


M 


{ fir)- lOAf (cm) 
, m 3 ) A(JQt?) 


(17) 


where M is the cloud liquid and 
A is the extent of the cloud. 
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Figure 2. Surface Measurements of Atmospheric Parameters During 
the Cloud Event . 


The GSW predictions are a 
little closer to the Slobin 
models but also exhibit some 
disagreement during the first 
half of the time period. 

All of this probably 

points to some physical 

phenomena that is not being 
accounted for in these simple 
"state equation" models. 

Improvements in modeling the 
vertical temperature profile, 
for example, might help 

matters. We are currently 
using simultaneous measurements 
of vertical temperature 
gradients and cloud liquid to 
improve the model . 

It is also of interest to 
note that the GSW model 
predicts the lifting 
condensation level to be a 

kilometer or so below the zero 
degree isotherm as shown in 


figure 6. We are now analyzing 


measurements of the lifting 
condensation level to improve 
cloud base altitude 
predictions . 

A striking result is shown 
in figure 7. Although the 
Altshuler and Dintelmann- 
Ortgies models were derived 
quite differently, they predict 
almost identical cloud 
attenuation time series 
patterns during the cloud 
event. Note however that the 
absolute magnitudes and the 
dynamic range of the patterns 
do differ. 


8 . Continuing Work 

The complexity of cloud 
physics and the lack of 
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Figure 3. Cloud Liquid and Vertical Cloud Extent Predicted by 

Dintelmann-Ortgies Model and NOAA Measurements of Cloud 
Liquid. 




Figure 4. Cloud Liquid and Vertical Cloud Extent Predicted by GSW 
Model and NOAA Measurements of Cloud Liquid. 
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Figure 5. Slobin Cloud Models [Slobin, 1982] . 


Vertical Temperature Profile Vertical Water Vapor Density Profile 



Temperature (deg C) Surf.Temp.— 22-06 C ; Surf.Density— 12.83 g/m A 3 

Figure 6. Vertical Profiles of Temperature, Saturation Water Vapor 
Density, and Actual Water Vapor Density as Predicted by 
the GSW Model. 
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Dintelmann-Orlgies Altshuler -Marr 



Figure 7. Comparison of Zenith Attenuation Predicted by the 
Dintelmann-Ortgies and Altshuler-Marr Models. 


measured data has always 
hampered cloud liquid research. 
Now as data begins to trickle 
in, we are seeing the 
beginnings of a new cloud 
liquid science — a blend of 
theory and experiment . The 
models presented here are a 
building block toward the 
understanding of cloud 
attenuation. As we continue 
working with more data sets at 
various locations, we are 
seeking to improve temperature 
profiling and condensation 
level predictions. Gradually 
we hope to incorporate and 
validate more detailed cloud 
physics to describe the 
condensation and mixing 
proces se s associated with 
clouds. We openly welcome your 
critiques and ideas. 
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